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I. INTRODUCTION 

Scattering of charged particles is an important and difficult topic in quantum mechanics. 
Even the two-body case is exceptional since it does not fit into the conventional time- 
dependent scattering theory approach. Several methods are proposed to overcome the diffi- 
culties caused by the long range behavior of the Coulomb interaction. In the time-dependent 
framework the formalism of Dollard^ modifies the Moller wave operator. In stationary ap- 
proach the method of van Haaringen^ introduced Coulomb asymptotic states and Mulherin 
and Zinnes^ promoted distorted asymptotic states. The connection of the latter two formal- 
ism was also studied^. Following references'"^ we will call the asymptotic state of Mulherin 
and Zinnes as Coulomb-distorted plane wave (CDPW). 

The recently developed surface integral formalism of quantum scattering theory'"- is 
valid not only for two-body scattering but it can be extended to scattering of three charged 
particles. Another advantage of this formalism is that it handles the short and longe range 
interactions on an equal footing. It is said in' "All the results ... rely on the asymptotic forms 
of the plane wave and the Coulomb-distorted plane waves" . The CDPW indeed plays an 
important role in the surface integral formalism. In this framework the full scattering wave 
function is composed of as a sum of the CDPW and a so called scattered part. Furthermore 
the surface integral expression of the scattering amplitude also refers to the CDPW. 

A three-dimensional leading order asymptotic form of the CDPW was given in'. Mathe- 
matical interpretation of this result was presented in reference' in section 5. The asymptotic 
expansion can be interpreted in distributional sense and the properties of the considered test 
function space plays an important role. Based on our convergent partial wave (pw) expan- 
sion of the CDPW and the asymptotic form of the pw components we will derive a three- 
dimensional leading order asymptotic form for the CDPW. Our expression is in agreement 
with the result of Ref.-'' if the test function space given by Taylor- is used. 

Interestingly the pw expansion of the CDPW has been given only recently^. In this work 
the pw components are applied for the description of scattering process of two particles 
based on the complex scaling method. Here we further study the pw decomposition and 
derive four different analytical forms. The asymptotic expansion of the pw component given 
in^ contains a recursion relation for the expansion coefficients. In this paper we derive an 
explicit expression for them. 



Our results concerning the CDPW are presented at Chapter [TTl Since the generahzed 
hypergeometric function 2-^2(0, a;a + I + l,a — I; z) plays a substantial role in the formal- 
ism therefore in the Appendix the properties of this function are examined in detail. The 
summary is given at chapter IIIII 

II. PARTIAL WAVE EXPANSION AND ASYMPTOTIC FORMS 

We write the pw expansion of the Coulomb-distorted plane wave in the form 

00 
e'^'^ikr T kr)±'^ = ^(2/ + l)r/^^(7, kr)Pi{cos{^)), (1) 

1=0 

where the upper and lower signs correspond to the post and prior form of the CDPW. The 
angle between the vectors k and r is signed by •& and the Legendre polynomial is denoted by 
Pi{x). The real variable 7 in ([1]) is the Sommerfeld parameter and k is the momentum vector. 
The bold face letters k and r denote elements of the three-dimensional Euclidean space M^, 
while k and r denote the corresponding magnitudes. The radial function r^ (7, kr) can be 
calculated by the integral 

r/^)(7, kr) = -{kr)^'^ / e±^^'^^(l T x)^'"'Pi{x)dx. (2) 

A compact expression for the pw component r^ (7, kr) can be given for arbitrary /. Using 
([2]) and the formula 2.17.5.6 in-- we get 

r/+)(7, kr) = . ^~'^} {2krY^'e'^' 2^2(1 + 27, 1 + ^7; / + 2 + 27, 1 + ^7 - /; -2ikr) (3) 
(l + z7ji+i 

and 

r/-)(7,fcr) = (-1)' ^ ^'^\ (2A;r)-^^e-^^%F2(l-z7, 1-^7; / + 2-27, 1-27-/; 2zA;r), (4) 
(1 -11)1+1 

where (a)„ is the Pochhammer symbol and 2-^2 is the generalized hypergeometric function^i. 
From the the explicit forms ([3]) and (jlj) one can notice that r^ (7, kr) = (— 1)V^ (7, kr)*. 
In order to further study the function r^ (7, kr) the properties of the function 2-^2(0-, a;a + 
I + l,a — I] z) are investigated in the Appendix. In the case of the post form of CDPW the 
following identifications can be made a = 1 + 27 and z = —2ikr and in the prior form we 
have a = 1 — i'-f, z = 2ikr. Using Theorem 1 of the Appendix and the formula ([3]) we can 



write down three equivalent expressions for the partial wave component 

^r^(7, kr) = ^^^^ EU(-1)"0^ 7(1 + n + n, 2^kr){2^kr)-- , (5) 

r/+)(7, kr) = {-iy{2kry^e^''^ ELo(-I)" (ifS^ 1^1(1 + H, 2 + ^7 + n; -2tkr) (6) 

and 

r/+^(7, A;r) = ^ (e*'=''/t+(l + ^7, -2ikr) + e-*'^'^/t;"(l + 27, -2ikr)) . (7) 

In charge-less case (i.e. 7 = 0) the CDPW goes into an ordinary plane wave. Us- 
ing Propositions 2 of the Appendix we get the expected result for the pw components i.e. 
Tj (0, kr) = i^ji{kr) i.e. the post and prior forms are identical. 

Using Theorem 2 of the Appendix we get for the asymptotic expansion of r^ (7, kr) as 

r — )• 00 

pikr 

r/+)(7, kr) -i^e^^'^ r(l + i^) ,F, (l + ^7, -/, / + 1; 1; {2tkr)-') 

^-ikr °o on (8) 

(2A;r)^-(-l)'5^4')- 



2ikr ' ' ' ^ Z^ n ^_2ikr)^ 



and 



„ -ikr 



r/-)(7, A;r) ~ - (-l)'^^e^-/2 r(l - i^) ,F, (l - Z7, "/, ^ + 1; 1; {-2ikr)~') 

pikr °^ r\n \^) 



+ k:(2*'-)""EK 



The expansion coefficients dn satisfy the recursion (lA.20tlA.2T]) and the solution of the 
recursion is given by (IA.22p . In these equations we have to identify a by 1 + 27 (1 — i'^) 
for the post (prior) form of the pw component. The expressions ([8]) and (|9]) show both the 
dominant and sub-dominant terms and the nature of the terms are determined by the sign 
of Im(iA;r). Instead of the study of the explicit form of r^ (7, kr) the asymptotic expansions 
dH]) and (jH]) can be obtained directly from ([2]) with a straightforward but tedious calculation 
based on the theorem 227 of Ref.— . 

Now we derive a three-dimensional leading order asymptotic form for the CDPW. Sub- 
stituting expressions ([8]) and (|9]) into ([1]) and keeping only the leading order terms we get 
our result in the form of a distribution 

(^^""{kr T kr)±^^ ~ ±-|^ (e^'^^e^^'^T{l ± i7)(5(f ^ k) - e^'^'' {2kr)^'^ 8{v ± k)) . (10) 



It is known that the prior form of the CDPW can be simply obtained from the post form 
expression: one has to replace k with — k and make complex conjugation. According to f fTOj) 
this procedure is valid also for the asymptotic forms. 

In charge-less case the three-dimensional asymptotic expression fITOl) goes into the well 
known form^^ 

e^kr _ ^ (e''''-6(r - k) - e-'^''6{r + k)) . (11) 

ikr \ / 

On the test function spaces D"^ of reference^ our result goes into equations (197) and (198) 
of^ 

(^^""{kr T kr)^^'^ ~ T^e^^^"(2fcr)±^^5(f ± k) on D^. (12) 

ikr 

On test functions belong to D^ the contribution of the first term of fITOl) is zero. Our 
rigorously derived result ( [T0|) differs from the expression of Ref.->^. However the surface 
integral formalism, as explained in reference, depends only on the validity of (IT2l) . 



III. SUMMARY 

The generalized hypergeometric function 2-^2(0, a;a + l + l,a — l; z) is expressed as a finite 
sum of confluent hypergeometric functions. Three equivalent expressions are determined 
and so the numerical evaluation can be simplified. Dominant and sub-dominant terms are 
identified in the asymptotic expansion. A closed expression is given for the coefficients 
of the asymptotic expansion. Using the previous results four equivalent expressions are 
determined for the partial wave component of the Coulomb-distorted plane wave. The 
asymptotic expansion of the partial wave component is deduced. Rigorously derived leading 
order three-dimensional asymptotic form of the Coulomb-distorted plane wave is given. 
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Appendix: The hypergeometric function 2^2(0, a;a + I + l,a — I; z) 

We will express the function 2-^2(0, a;a + l + l,a — l;z) in terms of finite linear combination 
of confiuent hypergeometric functions. In order to do this we need a proposition. 
1. Proposition. Let / be a non-negative integer and let X; := {x G Z | x < /}. If 
z e C\{0} and a G C\Xi then 

^ / / \ —k 

$^(-1)' ( ^ ) (/+l)fc(l-a)fc iFi(-fc, a-A;; z)^ = (-1)' 3Fi(l-a, -/, /+!; 1; -1/z). (A.l) 

fc=0 ^ ^ 

Proof. The definition of the confiuent hypergeometric function can be turned into the form 

k 






,=0 .sJia~k), 



Using this expression we get for the the left hand side of flA.l|) 

/ k 



EE(-i: 



k+sl ^\ fk\ {I + l)fc(l - a)fc^,_fc 



,kj \s J k\(a — k) 

k=0 s=0 \ / \ / \ I 

We can rearrange the summation indexes and get 

I 



s=0 

where 



Q— n ^ / 



' ' kfi-A ii + i)i-ks\ 



. k ){l + l)s{l-k)V 

Using the formula 0.160.2 of-- we get Si{s) = 1 and we can recognize that the summation in 
fTOj) gives the function 3Fi(l - a, -I, / + !;!; -l/z). D 

1. Theorem. Let / be a non-negative integer and let z G C\{0}, a G C\Z then the function 
2-^2(0, a]a + l + l,a — l;z) can be given in the following forms 

(-i)'g^ EU(-i)1D ggr ^Fi(«' <^ + k + i;z), (A.3) 

i-^)-^^ ELo (D ^ 7(a + k, -z) {z)-^ , (A4) 

{^(-z)- (/.+ (a,^) +e^Kr(a,^)) . (A.5) 

Here we introduced the notations 



n\{l — n)\ z^' 



KTia, z) = (-1)'+^ Y. 'ui-nV J U{-a + l,-a-n + l; -z) (A.6) 



n=0 



and 



/€+(a, z) = T{a) sFi{a, -/, l + l;l; -1/z). (A.7) 

The confluent hypergeometric function of tlie second kindil and the incomplete gamma 
function^, are signed by U{a, b, z) and 7(0, z), respectively. 

Proof. First we show that flA.SP equals flA.4p . Using the identity 5.3.5.4 of^'' and the 
expression 7.11.3.1 of-° we get for (1A.3I) 

I k 



-Z 



,-a Wl+l 



/' h—n Q=n \ / \ 



l\fk\{l + l] 



k=0 s=0 

Rearranging the summation we have 

I 



k\ 



7(a + s, -z)(-2;)' 



-^y 



a Wm ^(_i).(^^^ ([+11, ^(, + ,^ _,)(_,)-. ,^(,) 



1 - a 



s=0 



s\ 



and using that ei{s) = 1 we get ( 1A.4I) . 



Next we show that 2-^2(0, a; a + / + 1, a — /; 2;) equals (1A.4I) . Using the 



p 



(ap_i), a + I 
(bq-i), a 






(A.8) 



identity, which was proved inJ^, we get for 2-^2 («, a; a + / + 1, a — /; 2) 



A:=0 s=0 ^ 



l\ fl\ {a)i+i 



kj \sj l\(a — /); 



7(a + /c + s, -z){-zy 



(A.9) 



The following expression, valid for k > 0, can be easily established using the formula 8.356.9 



o£i^ 



fe-i 



7(6 + A;, ^) = {b)kl{b, z) - e-'z'' J^ 3^z^ 



(A.IO) 



and for /c = we have 7(6 + 0, 2;) = {b)o'y{b, z) = 7(6, z). Choosing 6 = a + s in ( JA.IOI) and 
substituting it into (1A.9P we get 



(-^) 



-a («) 



l+l 



n 



/\ r(a + k)T{l + m + l- k){-z) 



I k-l 



EE(-i) 



1 



(A.ll) 



kJ {a-l)kT{l + a + m + l) r{-[k - 1 - m])' 



k=l m=0 

In the derivation of this expression we used 0.160.2 of^^. We have to consider the case 
a G C\Z and I > k > 1 so we have that k — 1 — m is a. non-negative integer. We get that 



the double sum in (lA.lip is zero and this means that 2-^2(0, a;a + I + l,a — I; z) is identical 
with f lXij) . 

Finally we show that flA.Sp equals flA.Sp . Applying the expression 7.2.2.2 o&^ we get for 



dAU 
where 



iFi(a,a + l + n;2) + (-2)-"r(a)Q(a,z) , (A.12) 



ei{a, z) =3-Fi(a, -/, / + 1; 1; -1/^) 



/ 7 \ — n 

According to Proposition 1 e;(a, z) = and so we derived flA.SP from flA.5l) . D 

In order to show that the pw component of CDPW is proportional to the spherical Bessel 
function in the case of charge- less particles we will need the following statement. 
2. Proposition. Let / be a non-negative integer and let z G C\{0} and a G C\Z then 

lim (1 - a)i 2^2(0, a; a + / + 1, a - /; ;z) = i\l + 1)! e''''^ji{iz/2). (A. 13) 



a— >1 



Proof. It is easy to establish that (l)/+ie^fi:^"(l, 2)/^ = -«'(/ + l)\e''/'^hf\iz/2)/2 and 
{l)i+iKJ{l,z)/z = -i\l + l)\e'''^hf\iz/2)/2 using 10.1.16 and 10.1.17 ofii. From these 
facts and from flA.SP the claim follows. D 



2. Theorem. Let / > integer, a G C, < Re(a) < / + 2 then the asymptotic expansion 
of the function 2-^2(0, a; a + / + 1, a — /; z) at \z\ — )■ 00 can be written in the form 

2-r2(o, a; a + f + 1, a — /; z - 




.^N ,. - n „^ .„^ / _/ _; _. ^ . (^-14) 



n=0 



Proof. We will use the expression (lA.Sp to determine the asymptotic expansion. First we 
notice that K^(a,z) is a finite sum, it is a polynomial of order / in the variable 1/z. This 
means that !'"_+| {—z)~"-Kl{a, z) can be considered as an asymptotic expansion at \z\ — > 00 
with respect to the asymptotic sequence^ z~°- j z^ . The function H^{a^z\ according to 
equation (IA.6p . is a finite linear combination of confluent hypergeometric functions of second 



kind and since the asymptotic expansion of U{a, b, z) is given by 15.5.2 of^^ we can properly 
rearrange the summation indexes and we can write 



/€; [a,z) ~ (-1) 



oo / min(n,i) 

E E (-1)" 



{l + m)\ (m + l)„_m(l - a)„_„ 



m\{l — m)\ 



[n — m)\ 



(A.15) 



n=0 \ m,=0 

Using the formulas 7.2.3.15 and 7.3.5 of^° we can carry out the summation over m on the 
right hand side of f lA.15P and get 

' -/, -/, 
1, 1 — a — I 



Ki [a, z) ~ 



^y (-^)° y^ (l-a)„(a); ^ 



-n 



' t^o («-^)' 



1 



(A.16) 



D 



The expressions ( 1A.14I) shows that if \Ta.{z) > than the dominant term is given by the 
second term of the right hand side of flA.14p and if lm(z) < then the roles of the terms in 
(]A.14p interchange. 

The complete asymptotic expansion of the function 2-^2(^1, ^2; &i, ^2; z) was given by Luke 
in chapter 5.11.3 of reference^^. However in our case the numerator parameters are equal to 
each others and we can not use this result. We rewrite the function 2-^2(0, a;a + l + l,a — l; z) 
in terms of Meijer's G function (see 5.11.1(1) in^^) 



2^2(0, a;a + I + l,a — I; z) 



'lY(a. 



1 — a 



"^G^ I - 



1 — a, 1 — a 



0, 



-a — L I + 1 — a 



. (A.17) 



For the asymptotic expansion of the Meijer's G function we can use formula 5.10(10) of^^ 
and if we apply 5.7.(13-15) and 5.9.2 of Ref.— we get 



2-^2(0, a;a + l + l,a — l;z) 



-iy(a) 



:i-a) 



-H2o[-z) 



Taking into account 5.11.1(18) and 5.11.3(4) in^^ we can write 

z °° -1 



(A. 18) 



(A.19) 



n=0 



7(0 



The coefficients dn satisfy the recurrence relation (see 5.11.3(6) in 



15^1 



A{n + 1)4+1 = 2(2^2 - n{2a -3)-P-l-a + 1)4'^ -n{n-l-a){n + l + l- a)4-i (A.20) 



(0 



with initial conditions (see 5.11.1(18-20) in^) 



d: 



(0 



1 



1, 4" = ^(l-/^-/-a). 



(A.21) 



The asymptotic expansion (1A.18|) is valid in the region —^ + S< aig{—z) < ^ — 6 where 
6>0. 

We will show that our asymptotic expansion ( JA.141) and formula ( JA.181) are identical in 
the considered z region. Among the two terms at the right hand side of ( 1A.14I) the second 
one is the dominant in the considered z region and the first one can be neglected. To rewrite 
f lA.18|) and f lA.lOp we have to solve the recurrence relation ( 1A.20|) . 

3. Theorem. Let / and n be non- negative integers and let a G C, < Re (a) < / + 2 then 
the solution of the recurrence relation (]A.20|) with initial conditions ( 1A.21|) is 



rf« 



a)n{a)i 



2" (a - n)i 



3-^2 



L, -a-l + l 



(A.22) 



Proof. We prove this theorem by induction. With direct calculation it is easy to check that 



the claim is satisfied for n = 0,1,2. Assuming that (JA.221) is true for n = s and n = s + 1 
the claim, after some trivial simplification, is reduced to the validity of the equation 



[s + 2)(s + 2-/ -0)3^2 




[2(s + 1)2 - (s + l)(2a -3)-f -I 

. -I 

[s + l){s + 2 + I - a)^F2 



l + l 



l + l 



U 



This equation however is known to be true^. 

Using Theorem 3 and equation ( 1A.19I) we see that (1A.18I1 and the second term on the 
right hand side of ( 1A.14I) are indeed identical. 
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